Life times and chirality of spin-waves in antiferromagnetic and ferromagnetic FeRh: 
time depedent density functional theory perspective 
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The study of the spin excitations in antiferromagnetic (AFM) and ferromagnetic (FM) phases of 
FeRh is reported. We demonstrate that although the Fe atomic moments are well defined there is 
a number of important phenomena absent in the Heisenberg description: Landau damping of spin 
waves, large Rh moments induced by the AFM magnons, the formation of the optical magnons 
terminated by Stoner excitations. We relate the properties of the spin-wave damping to the fea- 
tures of the Stoner continuum and compare the chirality of the spin excitations in AFM, FM and 
paramagnetic (PM) systems. 
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Antiferromagnets (AFM) form a rich class of magnet- 
ically ordered systems characterized by the coexistence 
of non-zero local magnetization and zero net magnetic 
moment. In the contemporary physics, the interests in 
these materials has been refuelled by experiments on 
the ultra-fast laser induced magnetization reorientatior^, 
exchange-bias effect^ and the recent discovery of the ex- 
otic pnictide family of superconductors.'^ The latter ma- 
terials are often characterized by the proximity of the 
itinerant AFM and high-temperature superconductivity 
and the understanding of the spin excitations in the par- 
ent AFM compounds is crucial.'* 

Here, we study the spin-flip excitations in FeRh. There 
are important reasons for focusing on this system. It is 
antiferromagnetic at low temperatures and experiences 
the transition to the ferromagnetic (FM) state at 370 K. 
This allows the comparison of the AFM and FM spin ex- 
citations in the same material. Interestingly, the AFM- 
FM transition in FeRh can be initiated on the femtosec- 
ond time scale by the laser irradiationi^ 

The microscopic nature of the phase transformation in 
FeRh remains the topic of coirtroversial debates. It is 
widely recognized that the crucial role in the stabiliza- 
tion of the FM phase is played by the properties of the 
Rh atoms whose spin moments increase from zero in the 
ground AFM state to 1 /^b in the FM phase. Ju et al^ 
suggested that the driving force of the AFM-FM transi- 
tion is the appearance of the Rh moments in the AFM 
phase resulting from the fluctuations of the Fe moments. 
On the other hand, Gu and Antropo\^ put forward a 
principally different model where the appearance of the 
Rh moments in the AFM is not regarded and the phase 
transition is driven by the difference in the temperature 
dependence of the spin-wave entropy in both magnetic 
phases. Recently, we have shown that the spectrum of 
spin excitations in FeRh is more complex than considered 
previously." In particular we demonstrated that strong 
spin polarizability of the Rh atoms is the consequence of 
the implicit spin-polarization of the Rh electron states in 
the AFM ground state. 

All previous studies of the magnetic excitations in 
FeRh were performed within adiabatic approaches that 



map the system oir a model Hamiltonian of interact- 
ing atomic momentsj^^— A serious disadvantage of the 
method is the neglect of the one-electron Stoner excita- 
tions that can lead to such important effects as damp- 
ing, and even disappearance, of the spin waves as well 
as to strong modification of the magnon energies<^ The 
spin waves and Stoirer excitations are incorporated oir 
an equal footing is the calculations of the transverse dy- 
namic susceptibility within the framework of the time 
dependent density functional theory (TDDFT) i^"— Such 
studies for AFM systems are very scarce"''-'^^ because of 
the demanding character of the underlying computations. 
Recently, we have developed and implemented an efficient 
computational scheme for high accuracy calculations of 
the transverse magiretic susceptibility of complex mag- 
nets on the basis of linear response (LR) TDDFT J^^i^ 
Here, we report the study of the spin-flip excitations in 
AFM and FM phases of FeRh. 

On the experimental side, there is only one old mea- 
surement of the spin wave properties of FeRh performed 
by means of inelastic neutron scattering for both AFM 
and FM phases. The experiment revealed spin waves 
well defined in the whole Brillouin zone (BZ), in contrast 
to the spin-wave disappearance observed for large wave 
vectors in elemental id ferromagnetic bcc Fei^ii^ii^ Ad- 
ditionally, the experimental error bars for the spin-wave 
eirergies are very small, and in the case of AFM, the pre- 
sented error bars do not increase with increasing magnon 
momentum, contrary to the FM phase where certain in- 
crease is seen. In both AFM and FM phases only one 
spin- wave branch was detected. 

The calculation of the dynamic spin susceptibility 
within LR-TDDFT involves two steps,^iiiii^ Using the 
ground state Kohn-Sham electronic structure as the in- 
put, the non interacting magnetic response, XKs(q5'^)! 
of the system is found. The susceptibility does not give 
the physical magnetic response of the system, since the 
induced magnetization, leading to the additional induced 
internal exchange-correlation magnetic field described by 
the exchange-correlation kernel i^Txcr^, is not taken into 
account. Since the induced magnetization contributes 
to the effective external field and is, at the same time. 
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induced by this field, the problem must be solved self- 
consistently. The physical susceptibility x(q, w) is found 
from the Dyson equation 



X = (I - XKsi^xc) XKS- 



(1) 



The imaginary part of the susceptibility gives the inten- 
sity of spin-flip excitations with given energy uj and crys- 
tal momentum q. The information on the single-electron 
spin-flip spectrum, i.e. Stoner excitations^ is contained 
in the Kohn-Sham response, xks- The excitations form 
Stoner continuum. The spin-waves form at frequencies 
corresponding to the vanishing eigenvalues of the matrix 
I — XKS-f'^xc'^ If the spin- wave frequency appears outside 
the Stoner continuum, it features an infinite life-time. On 
the contrary, the magnons with energies within the con- 
tinuum decay due to the hybridization with single par- 
ticle excitations. This attenuation mechanism. Landau 
damping, dominates in metals. Inside the continuum, de- 
pending on its intensity, the physical picture of magnons 
varies between well defined quasi-particles (with life time 
much longer that their inverse energy) and spin- wave dis- 
appearance. We analyze the susceptibility by evaluating 
the loss matrix, — x^)/(2i). Its eigenvalues describe the 
intensity of spin-flip excitations, while the corresponding 
eigenvectors give the shape of natural modes of the sys- 
tem, i.e. spin-waves,i^ 
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FIG. 1: Schematics of band structures and Stoner continua 
in FM, AFM and paramagnetic material. 

This consideration is of a very general nature^^ and is 
valid for both AFM and FM and even for paramagnets 
featuring long-living collective spin excitations (param- 
agnons). There are, however, important qualitative dif- 
ferences in the structure of the Stoner continuum in 
AFM, FM and paramagnetic (FM) metals. They are il- 
lustrated schematically in Fig. [T] for a simple case of two 
bands corresponding to two opposite spin projections. 
The characteristic feature of the FM is the absence of 
the Stoner transitions in the low-w, low-q region. This 
is a consequence of the exchange splitting of the spin-up 
and spin-down bands. Magnons with small momenta are 
not Landau damped. On the other hand, in a metallic 
AFM the Stoner continuum starts at q = and w = 
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FIG. 2: Dispersion relation of spin-waves in the (a) AFM and 
(b,c) FM phase of FeRh. The dots (•) represent energy of the 
magnon peak maximum, while the error bars stand for the full 
width at half maximum (FWHM) of the peak. Diamonds (♦) 
denote the experimental spin-wave energies.^- (d) Example of 
spin-flip spectra in FM FeRh in the energy range of optical 
magnons. 



which follows from the spin degeneracy of the electron 
states. The AFM structure leads to the decrease of the 
BZ volume compared to the FM and FM cases and the 
formation of the second electron band. The transitions 
between states of the first and second bands form a sec- 
ond area of the high energy Stoner continuum. Also in 
PM crystal the states are spin-degenerate resulting in the 
low energy Stoner transitions. There is however an im- 
portant difference in the character of spin degeneracy in 
the AFM and PM crystals. In the AFM, the spin degen- 
eracy is the consequence of the presence of two equiva- 
lent magnetic sublattices and the wave functions of the 
degenerated states are shifted in the space with respect 
to each other. On the other hand, in the PM crystal 
the wave functions of the degenerate states are identi- 
cal. Because of the larger overlap of the initial and final 
states the intensity of the Stoner transitions in the low q- 
omega region tends to be stronger in the PM crystal. For 
real multiple-band systems the structure of the Stoner 
continuum is more complex, but the qualitative features 
depicted in Fig. [T] remain adequate. 

Now, we turn to the discussion of the results of our cal- 
culations. The obtained dispersion and damping of the 
spin waves of AFM FeRh are in good agreement with ex- 
periment, cf. Fig.[2Ji. The damping is not only very small 
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FIG. 3: Stoner continuum of AFM FeRh for small momenta. 
By selecting the q + Qap Fourier component of the KS sus- 
ceptibility, Qap = (1, 1, l)27r/aRh, we focus on the transitions 
between spin-degenerate bands. The spin-wave energies are 
marked with •. 
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FIG. 4: Transverse magnetization induced on the second Fe 
lattice and Rh site, assuming the deviation of the first Fe 
lattice is normalized to one. 
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over the whole BZ but it also clearly decreases with in- 
creasing momentum above 80 meV. In Fig. [2^ we do not 
show the values of full width at half maximum (FWHM) 
for the wave vectors close to F, since a reliable numerical 
estimation in this region is difficult. Instead, we make 
some qualitative remarks. The velocity of the electronic 
states at the Fermi level exceeds by an order of magnitude 
the magnon velocity given by the slope of the magnon dis- 
persion at the F point. This means that already for small 
momenta the magnons form inside the Stoner continuum 
and therefore for any q 7^ their damping is non-zero. 
The analysis of the Kohn-Sham susceptibility (Fig. [3]) 
shows that for small q the spectral density of the Stoner 
continuum has a shape of a narrow peak with its width 
vanishing in the 5 — > limit. This is a consequence of the 
decrease of energy differences between electron states in- 
volved in the transitions at small momenta. On the other 
hand, for increasing q the width of the peak increases 
and the position of the peak moves to larger frequencies. 
Since the energy width of the peak of the enhanced sus- 
ceptibility determining the position and life time of the 
spin wave cannot be larger than the energy area of the 
non-zero Stoner continuum the damping tends to zero in 
the g — >■ limit. Therefore, the spin-wave branch begins 
at g = with infinitely-long living magnons of vanishing 
energy. For q > the spin-wave peaks acquire non-zero 
width. On the basis of Fig. [3] we can conclude that the 
damping in the g-interval from to 0.1 first increases 
and than decrease again because of strongly decreasing 
spectral weight of the Stoner transitions at the spin wave 



FIG. 5: Chirality of spin-waves in different magnetic phases. 
Arrows denote the direction of magnetization precession. 

As mentioned in the introduction, the spin moments of 
the Rh atoms play very important role in the magnetism 
of FeRh. The calculation of the dynamic spin susceptibil- 
ity allows to investigate the spin-polarization of the Rh 
atoms by the AFM magnons. This question is in the di- 
rect relevance to the model of femtosecond generation of 
magnetization suggested Ju et al. ^ The analysis of the 
loss matrix shows that indeed such a polarization takes 
place. To understand this effect it is essential to recall 
that an AFM magnon breaks the equivalence of the two 
magnetic sublattices.— In order for the spins on both sub- 
lattices to precess coherently, the moments of one of the 
sublattices must deviate from the ground state direction 
stronger than the moments of the other. The ratio of the 
deviation angles of the two sublattices is q dependent. 
This is demonstrated in Fig. 2] where we show the trans- 
verse magnetization induced on the second Fe lattice and 
Rh site, assuming the deviation of the first Fe lattice is 
normalized to one. For g ^ the ratio tends to 1 what is 
a necessary condition for the fulfillment of the Goldstone 
theorem. On the other hand, at the BZ boundary the 
spins of the atoms of only one sublattice deviate from 
the ground state direction. The equivalence of the sub- 
lattices is restored by the presence of two degenerate spin 
wave branches corresponding to similar excitations with 
respect to different sublattices. 

In the AFM spin-wave the compensation of the Fe ex- 
change fields at the Rh sites is disturbed. This leads to 
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the appearance of non-zero atomic moments at the Rh 
sites. The value of the moments induced by the spin-wave 
is also q dependent, cf. Fig. E) It is zero for q ^> but 
assumes sizable values of 18% of the transverse Fe mo- 
ments for larger momenta. Although Fe atomic moments 
in FeRh are well defined, the richness of the physics of 
FeRh transcends the Heisenberg model incapable to de- 
scribe the magnetization induced on Rh sites. 

In the FM (Fig. ^ we have very good agreement with 
experiment in the low q region whereas closer to the zone 
boundary the deviation between theory and experiment 
becomes bigger. The damping has a clear tendency to 
increase with increasing momentum. A sizable damping 
is obtained for only rather large magnon wave vectors. 
This correlates with the property of the electronic struc- 
ture of the FM FeRh that has very small number of the Fe 
spin-up states in the nearest vicinity of the Fermi level.^ 

Since in FM the Rh atoms have a large moment of 
1 /iB, it is natural within the Heisenberg model to con- 
sider them as Heisenberg variables that leads to two spin 
wave modes: acoustic and opticaL-^ At higher frequencies 
we found sharp features in the enhanced spectral density, 
cf. Fig. [^Jl. The analysis of the eigenvectors of the loss 
matrix shows that the polar angle of the Rh magnetiza- 
tion deviation increases by 180° with respect to the di- 
rection of the Rh moment in the acoustic magnon. This 
is a characteristic feature of an optical magnon. The 
spectral feature, however, is characterized by a strong 
damping and appears at the energy position of abruptly 
increasing Stoner continuum, below the expected energy 
of the optical magnon. The resulting peak can be inter- 
preted neither as spin-wave nor as single-particle Stoner 
excitations, having a complex mixed character. The po- 
sition and width of the resonance as a function of wave 
vector is presented in Fig. [5] In principle, these res- 
onances should be observable in the inelastic neutron 
scattering experiment for sufficiently large neutron en- 
ergies and momenta. The peaks, however, carry much 
smaller spectral weight than the corresponding acoustic 
resonances. 

For comparison, we performed calculations for para- 
magnetic Rh. As expected the Stoner continuum starts 
in this case from the zero energy (not shown), similar 



to the case of AFM. It is however much more intense, 
in particular, because of the the stronger overlap of the 
wave functions of the states involved in the transitions. 
The enhanced susceptibility repeats the main features of 
the KS one. The enhancement in this material is not 
sufficiently strong to form a band of well defined para- 
magnon excitations similar to that found , e.g., in Pd.^ 
Since in paramagnetic metals the spin density vanishes, 
the symmetry with respect to the spin rotation is not 
spontaneously broken and the Goldstone mode does not 
form. In contrast to the spin waves in the AFM, the col- 
lective excitations in the FM are of longitudinal charac- 
ter and are connected with the formation of the non-zero 
magnetization instead of its rotation. 

The comparison of the spectral densities of the spin- 
flip excitations for a given q value in the cases of FM, 
AFM and FM reveals interesting similarities and differ- 
ences between them. In the FM, there is an asymmetry 
with respect to the change of the sing of the frequency: 
the spin-wave can be excited only at positive frequen- 
cies, which is a consequence of the broken time reversal 
symmetry of the FM ground state. In contrast, in the 
AFM the spectral density is symmetric with respect to 
the change of the sign of the frequency. These are two 
magnon branches that can be transformed to each other 
by the symmetry operation consisting of the product of 
the time reversal and space translation transforming one 
sublattice into other. By this transformation the sublat- 
tices interchange their roles. The paramagnet is invariant 
with respect to the time reversal that leads to the chi- 
rality property of the transverse magnetization similar to 
the case of AFM. 

To summarize, we show that although the Fe atomic 
moments are well defined in FeRh there are a number 
of important phenomena that are missed by the Heisen- 
berg model: Landau damping of spin waves, large Rh 
moments induced by the AFM magnons, the termina- 
tion of the formation of the optical magnons by Stoner 
excitations. AFM, FM and PM systems differ strongly 
with respect to the chirality of their collective spin-flip 
excitations. We hope this theoretical work will stimulate 
experiments on the spin-flip excitations in FeRh, e.g. us- 
ing modern high-flux neutron sources. 
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